Surface waves in solar granulation observed with Sunrise 
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ABSTRACT 

Solar oscillations are expected to be excited by turbulent flows in the inter- 
granular lanes near the solar surface. Time series recorded by the IMaX instru- 
ment aboard the Sunrise observatory reveal solar oscillations at high resolution, 
which allow studying the properties of oscillations with short wavelengths. We 
analyze two times series with synchronous recordings of Doppler velocity and 
continuum intensity images with durations of 32min and 23min, resp., recorded 
close to the disk center of the Sun to study the propagation and excitation of so- 
lar acoustic oscillations. In the Doppler velocity data, both the standing acoustic 
waves and the short-lived, high-degree running waves are visible. The standing 
waves are visible as temporary enhancements of the amplitudes of the large- 
scale velocity field due to the stochastic superposition of the acoustic waves. We 
focus on the high-degree small-scale waves by suitable filtering in the Fourier 
domain. Investigating the propagation and excitation of /- and pi-modes with 
wave numbers k > 1.41 /Mm we find that also exploding granules contribute to 
the excitation of solar p-modes in addition to the contribution of intergranular 
lanes. 
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Subject headings: waves — Sun: granulation — Sun: helioseismology — Sun: 
photosphere 



1. Introduction 



Solar p-modes are believed to be stochastically excited in the near surface layers of the 
Sun by turbulent convect ion. The pr ocess can be described as acoustic radiation by tur- 
bulent multipole sources (jUnno Ill964l ). Further theoretical descriptions of the excitation of 
sol ar oscillations follow t he wo rk on stochastic driving by turb ulent convection presented 
by iGoldreich fc Keeley I ( 119771 ) and iGoldreich fc Kumar I (119881 ). These models were suc- 
cessful in explaining the acoustic power distribution over a wide range in frequencies as is 
manifested in the ensemble of solar p-modes. 



Numerical simulations of convection, as described e.g. by lNordlund I ( 119851 ) or lStein fc Nordlund 

(120011 ). predict supersonic downdrafts which could play a role in the excitation of the 
solar oscillations. Ac t ually, the turbulent velocity field f rom t he numerical simulations 



by IStein fc Nordlund I (]200ll ) were used by ISamadi et al. I (120031 ) to determine the energy 



supply rate of convection into the solar oscillations. 

For the seismology of the Sun and sunlike stars, knowledge of the stochastic properties of 
the source function that drives p-modes would be a valuable input. In local helioseismology, 
e.g., the sensitivity functions which connect the acoustic properties of the p-modes with 
physical quantities in t he Sun depend on t he source function whose properties are not known 
and must be estimated (IBirch et al. 112004 ) . In addition, the Sun could serve as a paradigm for 
stellar seismology since the interaction of pulsations and convection and the energy gain and 



loss o f p-modes from radiation and convection is a general problem of stellar physics (IHoudek 
20061 ). 



When searching for the possible locations of the wave sources on the Sun, the turbulent 
downdrafts in the dark intergranular lanes are primary targets. So-called "acoustic events" , 
which are individual sunquakes with epicenters near the solar surface and located in the 
intergranular lanes, are assumed to feed continuously energy into the resonant p - modes 



of t he Sun. Studies of such acoustic events were carried out by iRimmele et al. I (119951 ) 



and lGoode et al. I (119981 ) with the result that such events have typical durations of 5 minutes 
and are caused by collapsing intergranular lanes. 



In this contribution our basis for further studies of the driving of the solar oscillations 



-3 - 



are Dopplergrams of the velocity above the s olar surface. Time series recorded with the 



IMaX instrument ( Martinez Pillet et al 



RISE (ISolanki et al. 



2010; 



Barthol et al. 



2010) aboard the balloon-borne observatory Sun- 



20101 ) provide the spatial and temporal resolution 



needed to measure the propagation of acoustic oscillations in the granules and the intergran- 
ular lanes. We assume that acoustic events could be observed as they eventually cause waves 
to propagate away from the location of the source. The unprecedented spatial resolution of 
the data makes it worthwhile to study the excitation and propagation of acoustic waves 
visible in the solar photosphere. As a new result we find that the acoustic power is not only 
enhanced in the intergranular lanes and along the borders of the granules but also above 
large granules, where splitting processes are about to occur. Surprisingly, the excitation of 
waves of high wave number seem to be a diagnostic of pending granule splitting. 



Observation and data reduction 



For our study we use a time series of solar granulation with a field of view of 51 ".5 
by 51". 5 located close to disk center. The data were obt ained with the Imaging Magneto- 
graph experiment - IMa X - ([M artin ez Pillet et al. 1120101) onboa rd the SUNRISE balloon- 



borne solar observatory (ISolanki et al. Il2010t iBarthol et al. Il2010l ) between 00:35:49 UT and 
00:59:57 UT (1st series) and 01:30:41 UT and 02:02:48 UT on June 9, 2009 (2nd series). The 
IMaX in strument and the reduction process which was applied to the data is described in 



detail in iMartmez Pillet et al. I (120101 ) and shall be briefly summarized in the following. 



The IMaX filtergraph was operated in the V5-6 observing mode, recording all Stokes 
parameters at A = [±80, ± 40, + 227] mA with respect to the line center of Fe I 5250.2 A 
(<7eff = 3, x — 0.121 eV). The exposure time for a picture at a single wavelength position 
was 6 s, while the entire line was scanned in 33 s. 

In the data reduction process, the dark current was removed, and a flat field procedure 
damped residual impurities, i.e. dust particles or scratches. In addition, the interference 
fringes - caused by the polarizing beam splitter - were removed from the images by atten- 
uating their contribution in Fourier space. Furthermore, a spatially dependent blue-shift 
over the FOV was calibrated and removed. Phase diversity (PD) measurements that were 
recorded close in time to the actual data, were used to reconstruct all pictures. The cor- 
rected images of the two cameras were merged to increase the S/N ratio which is of the 
order of 4 x 10 2 related to the continuum intensity (I c )in the reconstructed data. In a last 
step, crosstalk induced by residual motion was removed. Thus, we obtained time series with 
lengths of 24 and 32 minutes, resp., and a cadence of 33 s. A spatial resoluti on between 0."15 
- 0."18 was achieved by combining the on-board image stabilization system (IGandorfer et al. 
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20 id Berkefeld et al. Il2010l ; Martinez Pillet et al. 1120101 ) 



Besides Doppler velocity (v Dopp ) and continuum intensity (J c ), which we will use for 
further investigation, other quantities like line center intensity {I center), linear and circular 
polarization, and line width (FWHM) were computed. The quantities vd opp , ^center, and 
FWHM are inferred from a Gaussian fit to the four points in the line after they have been 
normalized using the continuum intensity at +227 niA. The velocity scale has been calibrated 
assuming a convective blue shift of 0.2 kms" 1 at disk center, and the absolute error in vo 
was estimated to be of the order of ±0.2 kms -1 . 



opp 



3. Solar Oscillations in IMaX Data 

3.1. Acoustic Power 

The two time series of the Doppler velocity measured by IMaX reveal a "flying shadow" 
pattern with a characteristic time scale of five minutes and a characteristic spatial length 
scale between 5" and 10". Su ch a movie will be par t of the electronic supplemenatary material 



of the introductory letter by lSolanki et al. I (120101 ) . This pattern is due to the solar /- and p- 



modes present on the Sun. The "flying shadows" originate from the stochastic superposition 
of a multitude of solar acoustic eigenmodes of various spatial wavelength, temporal frequency, 
and amplitude, which cause these "Five-Minute Oscillations". For a quantitative analysis 
we investigate the power distribution in the diagnostic diagram, where the spectral power 
density is displayed as a function of wave number k and temporal frequency u, cf. Fig. [TJ The 
frequency resolution in the Fourier domain is 0.52 mHz. We find a significant enhancement of 
oscillatory power between 2 and 7 mHz in the /- and p-mode area. Ridges of solar acoustic 
oscillations are at best visible only in outlines at high wave numbers. The /-mode ridge 
seems to extend until a wave number of ~ 5 Mm -1 , where the ridge is cut by the acoustic 
cut-off frequency at « 5.3 mHz (see Fig. [T]). The velocity of the superposition of /-modes 
has an rms-amplitude of 147 m/s with a maximum velocity of 902 m/s. From the location 
of enhanced power in this diagram we conclude that besides the superposition of large-scale 
standing acoustic waves with low wave numbers also the small-scale waves with high wave 
numbers are excited. These small-scale waves are expected to have a shorter lifetime than the 
large-scale waves. In the following we will focus on the study of these small-scale oscillations. 
Before doing so, we would like to note that besides the sonic power we find also subsonic 
power below the /-mode. These subsonic velocities correspond to the convective granular 
motions we observe. 
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3.2. Filter 

We study the excitation and propagation of solar p-modes by first filtering the data. In 
order to remove trends in the time series a "first-difference" filter is applied. Further filtering 
for the /- and pi-modes is carried out in the Fourier domain. A fast Fourier transformation 
is applied to the time series. The result of this is multiplied by a filter that acts as a 
band-pass including the /-mode and the pi-mode ridges. Figure [1] outlines (gray lines) this 
three-dimensional filter as projection onto the I — v diagram. The filter is based on Hanning 
windows with increasing widths as a function of the absolute magnitude of the wavenumber 
k. Further filter properties are that the filter allows only those modes with wave numbers 
1.4 Mm" 1 < I < 5.7 Mm" 1 to pass, i.e. the filter removes the contributions of large-scale 
waves and subsonic granular motions. After filtering, the data is transformed back into the 
time domain. The resulting time series contains only the small-scale oscillatory motions. 



3.3. Excitation and Propagation of Running Waves 

The filtered data allows the excitation and propagation of small-scale acoustic waves to 
be seen. For a quick investigation, acoustic power maps are determined by averaging the 
square of the filtered velocity time series over time and at each pixel. Figure [2] displays 
this temporally averaged acoustic power of the oscillatory velocity field after filtering. In 
this map the granular network is visible, since the intergranular lanes and the borders of 
the granules exhibit strong acoustic power. However, there is also significant power present 
above the granules. This is visible in detail in the enlarged view shown in the bottom panel 
of Fig. |2j Inspecting the complete series of Doppler maps, the origin of this acoustic power 
enhancements on top of the granules are oscillations which are either excited in the inte- 
granular lanes propagating through the granules or waves which are excited inside splitting 
or exploding granules. 

A time sequence of such an exploding granule is shown in Figs. |3] and HI During the 
observation the displayed granule continues to grow within the first half of the sequence 
before it fragments in the second half. Already at the beginning of the recording (t=199s) 
parallel wave fronts are rippling on top of the granule. Surprisingly, these wave fronts are 
aligned with an intergranular lane that is forming during the sequence but finally splits the 
granule 265 s later (t=464s). Further splittings of the granule occur at later times which 
seem to excite additional wave fronts on a smaller scale. Such events happen in multitude 
in the r ecorded time series and not only at this particular location. Based on a wavelet 



analysis iBello Gonzalez et al. I (120101 ) study the distribution of acoustic power in space and 



time. They find also increased power above small fast-evolving granules, splitting granules, 
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dark dots and forming lanes. 

Another finding, which can be observed in all granules of Fig. [3] and Fig. HI is that the 
wave amplitude itself is attenuated in the granules. The centers of stable granules show only 
weak small-scale oscillatory motions. After approximately four wavelengths the high-degree 
oscillations are damped away. Therefore, waves excited in the intergranular lanes or in the 
splitting process of granules seem to lose energy while probing a granule. This might explain 
the comparatively high acoustic power along the borders of the granules. 

4. Conclusions 

We have investigated the properties of small-scale solar p-modes as recorded by the IMaX 
instrument during the flight of the balloon-borne Sunrise observatory in June 2009. We 
analyzed the Doppler velocity and continuum intensity observations to study the potential 
of such data sets for investigations of the excitation and propagation of high-degree acoustic 
waves on the Sun. Based on two sequences of data, we find that there is a strong relation 
between increased acoustic p-mode power with high wave number above granules and the 
onset of the splitting process that fragments granules. The splitting of granules triggers the 
excitation of waves several minutes before the actual splitting becomes visible in continuum 
radiation. A possible explanation for the excitation of such waves could be related to the 
formation of new intergranular lanes and connected turbulent motions which set off acoustic 
waves. 

Therefore, we conclude that not only the vigorous motions in the established integranu- 
lar downdrafts are sources for acoustic waves, but also the processes related to the formation 
of new intergranular lanes. In future studies, e.g. based on either data from a second Sunrise 
flight or from other instruments in space or on ground, we plan to investigate the excitation 
process of solar p-modes on the basis of spatially highly resolved data in greater detail. For a 
more quantitative analysis we aim to obtain new contiguous observations that cover a longer 
time period, which shall allow to resolve the ridges of the p-modes in the diagnostic diagram. 
Such longer time series are needed to measure e.g. the Green's function which is the respose 
to an "acoustic event" on the Sun. This would help determining the statistical properties of 
the solar oscillation's source function, which is needed as an input to sensitivity kernels as 
they are calculated, e.g., for local helioseismology. Furthermore, it remains an open task to 
study the damping, transmission, and refraction of solar p-modes in granules. Observations 
of such quality could be useful for such investigations and might allow obtaining an acoustic 
image of the formation and the structure of granules. 
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Fig. 1. — Spectral power density of the measured Doppler velocity as a function of wave 
number and frequency in arbitrary units. The color code reflects the mode energy: large 
amplitudes are coded dark, amplitudes compatible with zero are coded white. The black 
solid line indicates the theoretical position of the /-mode. The gray solid lines mark the 
borders of a Hanning band pass filter, which is applied to the data. 
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Fig. 2. — Acoustic power maps of the filtered data. White locations on the Sun 

with high acoustic power. Top: Average acoustic power of the whole field-of-view. Bottom: 
Average acoustic power of an 8". 25 x 8". 25 section. The white rectangle in the top panel 
marks the area which is shown enlarged in the bottom panel. In the bottom panel only 
values below a threshold of 1500 (m/s) 2 are shown to make the acoustic power above the 
granules visible. 
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Fig. 3. — Time sequence of the filtered velocity time series (left) and the synchronously 
recored continuum images of the same area (right). The snapshots are separated by 66s. 
The exploding granule is located in the center. The spatial extent of this area is 8". 3 x 8". 3. 
The displayed area corresponds to the map shown in the bottom panel of Fig. [2J The series 
is continued in Fig. |H 
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Fig. 4. — Continuation of the time sequence shown in Fig. [3j 



